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Accurate ab initio calculations have been carried out to study the valence electron 
removal energies and oscillator strengths of astrophysically important electromag- 
netic transitions of quadruply ionized vanadium, V^^. Many important electron 
correlations are considered to all-orders using the relativistic coupled-cluster the- 
ory. Calculated ionization potentials and fine structure splittings are compared 
with the experimental values, wherever available. To our knowledge, oscillator 
strengths of electric dipole transitions are predicted for the first time for most of 
the transitions. The transitions span in the range of ultraviolet, visible and near 
infrared regions and are important for astrophysical observations. 

Subject headings: atomic data — stars: circumstellar matter — ISM: atoms 



1. Introduction 



Vanadium (V) is very important constituent in the atmosphere of dwarf stars (IKirkpatrick et al 



19991 ) as it condenses into solid solution with Ti-bearing condensates (as observed in mete- 
orites) and not as pure vanadium oxide (V2O5) as commo nly assunied. A bundance of V in its 
ionic or molecular form is used to constrain temperature (jLoddersll2002l ) at the dwarf transi- 
tion. The abundance of V has been seen for 46 normal stars of spectral types from G8 to MO, 
and their accurate estimates depend on the precise e valuation of transition amplitudes. This 
abundance is very important in Sun (jBonsacklll959l ) and the observed maximum abundance 
among similar stars declines with surface temperature. The abundance excesses of this odd 
Z element by 1.2 dex in the hot peculi ar star, like 'star 3 cen A' , has very interesting and 
important applications in astrophysics (ICowley fc Wahlgernll2002l ). The transition between 
low-lying energy levels is also of importance to the identification of V emission lines in solar 
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spectra. Spectroscopic studies of this mult i- charged ion are important in numerous apphca- 
tions in laser and plasma physics apa rt from astrophysics, li ke, the operations of short-wave 
laser plasma, plasma diagnostics etc (jVainshtein et al.lll986l ). 



The dominant configuration to the ground state of quadruply ionized vanadium (V^^) 
is [Ar]3(i3/2 (^-03/2), with strong correlations from 3d orbitals. This is unlike the isoelectronic 
neutral potassium (K), whose ground state is dominated by [Ar]45i /2 ('^S-\ /-?) . Exp erimental 



ionization potentials (IP) of groun d and excited states of this ion (IDeurzeru 119771 ) are tab- 
ulated in the NIST web database (INIST website reference! ) . The comparison highlights the 
effective considerations of quantum many-body theories, especially the correlation contribu- 
tions applied to evaluate these energy eigenstates. There is no line strength estimation among 
different low lying states obtained from the change of outer most valence orbitals of V^~^ to 
our knowledge. For few transitions, preliminary semi-empirical results a re available wherein 
exper imental energies had been used to calculate oscillator strengths (IMalinovskaya et al. 
1992h . 



Here we have employed the relativistic coupled cluster method with single and double 
excitations (RCCSD). This is one of the most po werful highly correlated appro aches due to its 
all order behavior for the correlation operator (ILindgren and Morrison! 119851 ). The all order 
behavior comes from the generation of correlated RCCSD states usi ng exponentia l structure 
of the cluster operators over the Dirac- Fock (DF) referenc e states (ICoester!!l958! ). which is 
explicitly e xplained in our recent paper (jSahoo et al.!!2004l ). The contributions due to Breit 
interaction (!Breit!!l929l ) which is four orders smaller than the Coulomb interaction have been 
neglected. For a univalent atomic system a set of coupled equations for the cluster operators, 
T„,, may be obtained from the Bloch equations fo r core and valence electrons separately 



(ILindgren and Morrison! ! 19851 : !Debasish et al.!!l989! ). Because of this consideration, though 
the cluster operators are only of single and double excitations kind, their combinations can 
produce partial triple excitations which nominate the method as coupled cluster method with 
single, double and partial triple excitations (CCSD(T)). Contributions fro m these partia l 



excitations in the ionization potentials (IPs) for other systems are well kno wn (!Bartlett 



1989); 



this technique has been described extensively in one of our recent papers (iSahoo et al.!!2004j ). 
Therefore, it can take into account of different electron correlations, like core correlation, 
core polarization a nd pair correlation, effects exhaustively for specific type of excitations 
JSahoo et al.ll2004h . 



The oscillator strength for a transition from |\E'j) to |\I'/) is 



D 



fi\ 



where AEfi is the energy difference between final and initial states and gi = 2Jj -|- 1 is the 
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degeneracy factor of the initial state with angular momentum Jj. The electric dipole moment 
matrix element Dfi is defined as 



D 



fi 



(*/|d|^..) 



where, 



(2) 
(3) 

d= er is the electric dipole moment operator in the length (L) gauge form. |$)s' are the 
DF reference states. Sim ilarly, transition matrix elements can be obtained by replacing 
d in velocity gauge form (jLinlll978l : IW. R. Johnson et al. II2002I ). The connected parts of 
this expression will contribute and hence we only compute those parts in our dipole matrix 
element calculations. 

In the astrophysical literature, one uses weighted oscillator strength which is the product 
of degeneracy of the initial state and the oscillator strength and is symmetric with respect 
to initial and final states, i.e.. 



gf={2J, + l)f,f = -{2Jj + l)fj,. 



(4) 



2. Results and Discussion 

We have started with the generation of DF orbitals cons tructed as linear combinatio n 



of Gaussian type orbitals (GTOs), given by Chaudhuri et al. (IR.K. Chaudhuri et al.ll2000l ). 
The number of single particle orbital bases are chosen for different symmetries based on 
the number of occupied orbitals and requirement of the transitions of interest. All the core 
orbitals are considered active in our calculations to make correlations more exhaustive, es- 
pecially the core correlation part. Table I shows the excitation energies and fine structure 
splittings of many low lying states with single unfilled orbitals of V'^'^. For ionization poten- 
tials, we find very good a greement with the experimental data obtained from NIST database 



(INIST website reference! ) for all the cases. Leaving out the 3^5/2 state which has rather low 
binding energy, the average deviation from the experimental values are about 0.1%. The 
experimental uncertainty of the allowed transition wavelengths in the NIST values are very 
small, at the level of 2nd to 3rd decimal places, which are of the order 10 cm~^ in excitation 
energies. The fine-structure splittings are of the same order. Therefore, it is important to 
do higher precession experiments for fine-structure splittings. 



The oscillator strength for transitions between different energy states depend linearly 
on energy, and quadratically on the electric dipole matrix elements (Eq. (1)). This shows 
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Table 1: Ionization Potentials (IPs) and fine structure (FS) splittings (in cm~^) of V'^^ and 
their comparison with NIST values. The percentage of differences of our calculated IP results 
compared to NIST results are shown in parenthesis on the side of CCSD(T) IP results. 
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the necessity to calculate transition matrix elements appropriately. We present the elec- 
tric dipole transition amplitudes in both length and velocity gauges in Table II. The good 
agreement between these two forms from the calculations using the same initial and final 
wavefunctions highlights the fact that the numerical algorit hm is robust. A detailed d eriva- 
tion and interpretation of different gauges can be fou nd in (IBethe and Salpeter 1 119571 ) . and 
in few recent references (jW. R. Johnson et al. 1120021 ). The little differences in the values of 
these two forms for particular transition amplitude may be improved by considering nega- 
tive energy states in the calculation. In few transition amplitudes, there are considerable 
disagreement between the length- and velocity-forms. The det ail analysis might be possible 
using the approach suggested by (ISavukov and Johnson Il2000l ). In table III, we have given 
the ab initio oscillator strengths (f-value) for 'ground to excited' and 'excited to excited' 
transitions, which are also astrophysically important. Here we are pr esenting the oscillat or 
strength in the length gauge due to its comparative fast convergence (jPas &: Idresslll990l ). 



3. Conclusion 



Highly correlated relativistic coupled cluster theory has been employed to study the os- 
cillator strengths of the astrophysically important electric dipole transitions. To our knowl- 
edge these are the first calculations of oscillator strengths for most of the transitions pre- 
sented here. All the transitions are in the ultraviolet, visible or near infrared regions. The 
good agreement between transition amplitudes obtained from length and velocity gauge ex- 
pressions highlight the accuracies of our computational method. This work will motivate 
astronomers to observe these lines of V'^'^ to predict the abundances of these species in 
astronomical bodies and experimentalists to verify our results. 
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Table 2: Comparisons between length and velocity gauge transition amplitudes. 



Terms 





Length 


Velocity 


4pi/2 


1 


.74965 


1. 


,75033 


5pi/2 


0, 
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0. 
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-0. 


,07370 


6pi/2 


-0, 


.09081 


-0. 
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.81273 


-0. 


,97714 


4pi/2 


0, 
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3. 
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Length 


Velocity 
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Table 3: Transition wavelengths and oscillator strengths of V^~^. 
Terms 
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0. 


.3993 


5p3/2 


1929, 


.13 


1918, 


.22 


0. 
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.99 


1717, 


.14 


1. 


.4974 


6p3/2 


822, 


.17 


822, 
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.64 


1871, 


.18 


0. 


,2289 


4/7/2 


1811, 


.42 


1872, 


.01 


4. 


,5612 


55-7/2 


1490, 


.13 


1540, 


.61 


0. 


,2304 


55-9/2 


1490, 


.10 


1446, 


.24 


8. 


,5927 


65-7/2 


992, 


.33 


959, 


.49 


0. 


.0395 


6^9/2 


992, 


.33 


975, 


.19 


1. 


.3638 


4pi/2 


550, 


.71 


550, 


.24 


0. 


.0018 


4p3/2 


554, 


.57 


554, 


.17 


0. 


,0002 


5pi/2 


2741, 


.48 


2747, 


.64 


2. 


,9279 


5p3/2 


2780, 


.96 


2784, 


.77 


0. 


,5827 


6pi/2 


3643 


.91 


3651, 


.63 


0. 


,0969 


6p3/2 


3610, 


.28 


3654, 


.49 


0. 


,0193 


4/5/2 


2610, 


.86 


2656, 


.77 


0. 


.7421 


4p3/2 


554, 


.37 


553, 


.95 


0. 


.0020 


5p3/2 


2775, 


.82 


2779, 


.17 


5. 


.2505 


6p3/2 


3618, 


.99 


3664, 


.18 


0. 


.1735 


4/5/2 


2606, 


.33 


2452, 


.41 


0. 


,0529 


4/7/2 


2577, 


.89 


2450, 


.98 


1. 


,0516 


4pi/2 


438, 


.76 


438, 


.33 


0. 


,0008 


4p3/2 


441, 


.22 


440, 


.83 


0. 


,0002 


5pi/2 


1207, 


.68 


1207, 


.87 


0. 


,0334 


5p3/2 


1215, 


.28 


1214, 


.99 


0. 


.0057 


6pi/2 


5295, 


.58 


5260, 


.35 


0. 


.3127 


6p3/2 


5368, 


.23 


5254, 


.42 


0. 


.0618 


4/5/2 


1181, 


.63 


1148, 


.12 


0. 


,0455 


4p3/2 


441, 


.14 


440, 


.83 


0. 


,0021 


5p3/2 


1214, 


.73 


1214, 


.98 


0. 


,0530 


6p3/2 


5360, 


.42 


5254, 


.26 


0, 


,5571 


4/5/2 


1181, 


.26 


1148, 


.11 


0. 


,0032 


4/7/2 


1175, 


.38 


1147, 


.80 


0. 


.0648 



